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Abstract. A search for exotic + baryon via + — > p + Kg decay mode in the NOMAD f M iV data is 
reported. The special background generation procedure was developed. The proton identification criteria 
are tuned to maximize the sensitivity to the <9 + signal as a function of xf which allows to study the + 
production mechanism. We do not observe any evidence for the + state in the NOMAD data. We provide 
an upper limit on + production rate at 90% CL as 2.13 • 10 -3 per neutrino interaction. 

Key words, neutrino interactions, strange particles, exotic baryons, pentaquarks 

PACS. 13.15.+g Neutrino interactions - 13.60.Le Meson production - 13.87.Fh Fragmentation into 
hadrons - 14.40.Ev Other strange mesons 
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1 Introduction 

In the last three years an intense experimental activity has 
been carried out to search for exotic baryon states with 
charge and flavor requiring a minimal valence quark con- 
figuration of four quarks and one antiquark (such states 
are often referred to as "pentaquarks" ). Searches for exotic 
baryon states have a ~ 30 year history, but a theoretical 
paper by Diakonov, Petrov and Polyakov [1] has triggered 
practically all recent activity. 

The LEPS Collaboration was the first to report the ob- 
servation of the + (uudds) state with positive strangeness 
[113]. Then confirmations followed from DIANA (ITEP) g], 
CLAS [5115117115] . ELSA (SAPHIR) [5], old (anti) neutrino 
bubble chambers data (WA21, WA25, WA59, E180, E632) 
reanalyzed by ITEP physicists [TO], HERMES [IT], SVD 
(IHEP) [H], COSY-TOF [13], LHE (JINR) Q3], HEP 
ANL - HERA (ZEUS) [T5J[TS]. A narrow peak in the in- 
variant mass distributions of pKg or nK + pairs with a 
mass of ~ 1530 — 1540 MeV/c 2 and a width of less than 
25 MeV/c 2 was observed in all these experiments with 
significances of 4-8 tr's. Searches for narrow pentaquark 
states were then performed in almost every accelerator 
experiment in the world, providing evidence or hints for 
a variety of pentaquark candidates: + , 3jE~, ++ , and 
(9°. However, after this initial flurry of positive results, 
negative results, in particular from high statistics exper- 
iments, started to dominate the field. As an example of 
a negative search we quote the HERA-B experiment at 
DESY [T7J that observed neither the + resonance in 
the pKg invariant mass distribution nor the EE (1860) 
(another member of the antidccuplet of exotic baryons) 
decaying to EE~n~. Also the BES Collaboration [15] re- 
ported no 0(1540) signal in ip(2S) and J/tjj hadronic de- 
cays to KgpK~n and KgpK + n. The PHENIX experiment 
at RHIC [TI5] has seen no anti-pentaquark 0~ in the de- 
cay channel K~n. Also the BABAR [15] and the CDF [21] 
experiments have provided no evidence for + . Possible 
explanations for such a controversial experimental situa- 
tion could be ascribed to specific production mechanisms 
yielding pentaquarks only for specific initial state parti- 
cles. However, the CLAS experiment at Jefferson Lab has 
recently reported the results of a new analysis of photon- 
deuterium interactions with a statistics six times larger 
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than the earlier event sample which showed a positive re- 
sult. In this new analysis no + peak was seen [35]. A 
review of the experimental evidence for and against the 
existence of pentaquarks is presented in |23] . 

This article describes a search for the lightest member 
of the antidecuplet of exotic baryons, + , in the decay 
channel + — > p + Kg from a large sample of neutrino in- 
teractions recorded in the NOMAD experiment at CERN. 

The paper is organized as follows. In Sec. d we give a 
brief description of the NOMAD detector, and of the NO- 
MAD simulation program (MC). In Sec. [3] we present the 
event selection criteria and the tools for K s and proton 
identification. In Scc.[3]wc describe checks of the proton 
identification procedure and discuss the expected invari- 
ant mass resolution of pKg pairs. We describe in detail our 
procedure for determining the shape of the background 
distribution in Sec. 2J Based on the background determi- 
nation procedure and proton identification, we then de- 
velop a strategy for a "blind" analysis of the + signal 
by finding the proton identification criteria which maxi- 
mize the sensitivity to the expected signal. This approach 
is presented in Sec. [5J In this section we present also the 
method to estimate the signal significance and check the 
analysis chain using the observed decays A — > p + n~ and 
Kg — > 7r + 7r~. Finally, we "open the box", i.e. examine the 
signal in the data. The conclusions are drawn in Sec. [6] 



2 The NOMAD detector 

The large sample of neutrino interactions, about 1.5 mil- 
lions, measured in NOMAD together with the good re- 
construction quality of individual tracks, offer an excellent 
opportunity to search for 0+ -> p+K%. The NOMAD de- 
tector [21] consisted of an active target of 44 drift cham- 
bers, with a total fiducial mass of 2.7 tons, located in a 
0.4 Tesla dipole magnetic field, as shown in Fig. [1] 

The drift chambers |25j , made of low Z material (mainly 
Carbon) served the double role of a nearly isoscalar tar- 
get for neutrino interactions and of the tracking medium. 
The average density of the drift chamber volume was 0.1 
g/cm 3 . These chambers provided an overall efficiency for 
charged track reconstruction of better than 95% and a mo- 
mentum resolution of approximately 3.5% in the momen- 
tum range of interest (less than 10 GeV/c). Reconstructed 
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tracks were used to determine the event topology (the as- 
signment of tracks to vertices), to reconstruct the vertex 
position and the track parameters at each vertex and, fi- 
nally, to identify the vertex type (primary, secondary, V , 
etc.). A transition radiation detector [2"6"I|2"T] placed at the 
end of the active target was used for particle identification. 
A lead-glass electromagnetic calorimeter [28|f29] located 
downstream of the tracking region provided an energy res- 
olution of 3.2%/ y/E\GeV] 1% for electr omagnetic show- 
ers and was crucial to measure the total energy flow in 
neutrino interactions. In addition, an iron absorber and 
a set of muon chambers located after the electromagnetic 
calorimeter were used for muon identification, providing 
a muon detection efficiency of 97% for momenta greater 
than 5 GeV/c. 

Neutral strange particles were reconstructed and iden- 
tified with high efficiency and purity using the V°-like sig- 
nature of their decays [301131"] . Proton identification needed 
further development for the search presented here using 
information from the drift chambers, transition radiation 
detector and electromagnetic calorimeter. 

The NOMAD Monte Carlo simulation (MC) is based 
on LEPTO 6.1 [3J and JETSET 7.4 [33J generators for 
neutrino interactions, and on a GEANT [34] based pro- 
gram for the detector response. The relevant JETSET pa- 
rameters have been tuned in order to reproduce the yields 
of strange particles measured in CC interactions in NO- 
MAD [31] . To define the parton content of the nucleon for 
the cross-section calculation we have used the parton den- 
sity distributions parametrized in [35] . 

3 Event Selection 

We have analysed neutrino-nucleon interactions of both 
charged (CC) and neutral current (NC) types. These events 
are selected with the requirements: 

- The reconstructed primary vertex should be within 
a fiducial volume (FV) defined by \x,y\ < 120 cm, 
5 < z < 395 cm (see Fig. [T] for the definition of the 
NOMAD coordinate system) 

- There should be at least two charged tracks originating 
from the primary vertex; 

- The visible hadronic energy should be larger than 3 
GcV. 

The CC events are identified requiring in addition: 

- The presence of an identified muon from the primary 
vertex. 

The NC sample contains a contamination of about 30% 
from unidentified CC events. However, we do not apply 
further rejection against this background in order not to 
reduce the statistics. The event purity for the CC se- 
lection is 99.6%. The total sample amounts to about 1.5 
million neutrino events (sec Table [T]). 

3.1 Kg Identification 

K% mesons are identified through their V°-like decay 
K s — > 7r + 7r - using a kinematic constrained fit [30]I31|. With 





CC 


NC 


CC+NC 


N obs 


785232 


393539 


1178771 


N corr 


1017664 


481269 


1498933 



Table 1. Statistics of observed (N oos ) and efficiency corrected 
(N C orr) neutrino CC and NC events in the data. 



a purity of 97% we identify 15934 and 7657 K s mesons in 
the CC and NC samples respectively, thus yielding a to- 
tal statistics of more than 23k Kg's. The reconstructed 
Kg — > 7r + 7r~ invariant mass distribution in the CC 
(left) and NC (right) subsamples are shown in Fig. [2] 
The two distributions have the same K s mass mean value, 
497.9 MeV/c 2 , in agreement with the PDG value, and a 
width compatible with the expected experimental resolu- 
tion of ~9.5 MeV/c 2 . 



3.2 Proton Reconstruction 

The identification of protons is the most difficult part of 
the present analysis. As the NOMAD experiment does not 
include a dedicated detector for proton identification, we 
developed a special procedure for this purpose. The main 
background in the proton selection is the 7r + contamina- 
tion since pions are about 2.5 times more abundant than 
protons. However the 7r + contamination can be suppressed 
exploiting the differences in the behaviour of protons and 
pions propagating through the NOMAD detector. We use 
three sub-detectors which can provide substantial rejec- 
tion factors against pions: 

1. The Drift Chambers (DC). A low energy proton ranges 
out faster than a pion of the same momentum. Thus 
a correlation between the particle momentum and its 
path length can be used discriminator between 
protons and pions. The momentum interval of appli- 
cability of this method is below 600 MeV/c. 

2. The Transition Radiation Detector (TRD). The en- 
ergy deposition of protons and pions in the TRD is 
very different due to the larger proton ionization loss 
for momenta below 1 GeV/c, allowing a good pion- 
proton separation in this momentum interval. A mod- 
est discrimination is also possible for momenta above 
3 GeV/c because of relativistic rise effects. 

3. The Electromagnetic Calorimeter (ECAL). The proton 
sample can be cleaned further by taking into account 
the different Cherenkov light emission of protons and 
pions of the same momentum. 

As a preliminary quality selection for the candidate 
proton track we require : 

- more than 7 hits on the DC track in order to have a 
reliable fit; 

- the distance between the primary vertex and the first 
hit of the track to be smaller than 15 cm; 

- the relative error on the track momentum to be smaller 
than 0.3; 
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- the coordinate of the last hit of the stopped particle to 
be within a fiducial volume (FV) defined as \x\ < 120 
cm, -110 < y < 100 cm, 35 < z < 380 cm, which 
is smaller than the FV used in the event selection in 
order to reduce edge effects; 

- the particles reaching the TRD to cross at least 6 TRD 
planes (out of 9 in total) in order to allow for the TRD 
identification algorithm. 

The initial sample of positively charged tracks is split 
into two subsamples: those that passed the quality crite- 
ria and those that did not. The identification cuts can be 
applied only to the positive tracks that passed the qual- 
ity criteria. Charged particles were all assumed to be pi- 
ons in the standard reconstruction program, resulting in 
a systematic underestimation of the reconstructed pro- 
ton momentum (p rec ) at small momenta. The difference 
Ap between the true and reconstructed proton momentum 
was studied with the help of the MC, and the following 
parametrization was found to describe the effect : 

Ap = 0.33 • e "%^r GeV/c. (1) 

The difference Ap decreases with p rec , becoming neg- 
ligible at about 0.8 GeV/c. The reconstructed proton mo- 
mentum was corrected accordingly. The effect of this pro- 
ton correction was tested on a reconstructed sample of A 
hyperons identified by their displaced decay vertex. Fig. [3] 
displays the mean value of the invariant pir~ mass as a 
function of the reconstructed proton momentum without 
and with the Ap correction. There is an improvement 
in the reconstructed A mass when correcting the recon- 
structed proton momentum, especially at low momenta. 

4 The background 

Random _ftT^-proton pairs produce combinatorial back- 
ground in the Kg-p invariant mass (M(K s p) = M) dis- 
tribution. Understanding the shape of this background is 
crucial in the search for a possible + signal. We studied 
this background in three different ways: 

1. MC events contain no <9 + and could be used, there- 
fore, to study the background for this analysis. How- 
ever, the small fraction of proton-ifg pairs with an 
invariant mass in the interesting mass region would 
require a very large sample of MC events to reduce 
statistical fluctuations. 

2. We combined protons and K s s from different events 
in the data, thus making fake pairs, paying special at- 
tention that the original data distributions of multi- 
plicity, proton and kaon momenta, and their relative 
opening angle, were well reproduced in the final fake 
pair sample. 

3. A polynomial fit to the M distribution of the data 
themselves, excluding the c9 + mass region, can also be 
used to describe the background for the c9 + search. 

The fake pair technique is extensively used in the liter- 
ature. However, it is necessary to ensure that the two in- 
dependent events used in the mixing have similar hadronic 



jet momenta both in magnitude and direction. If two events 
with different jet momenta are mixed, then the fake pair 
technique systematically underestimates the background 
at small M(K s p) invariant masses. This is illustrated in 
the left panel of Fig. 2] which shows the Kgir + invari- 
ant mass distribution in the data, superimposed with the 
background generated without changing the hadronic jet 
directions. The background is normalized to the data at 
M(K%ir+) = 1400 MeV/c 2 . The data show a clear K*+ 
peak at ~ 890 McV/c 2 . However, the background under 
this peak is obviously underestimated. Therefore, in our 
procedure we first rotate each data event such that the 
hadronic jet momentum is aligned along the z-axis. We 
then select events respecting the original multiplicity of 
positive tracks and Kg in the data and make random 
pairs. The right panel of Fig. 2] shows the same K s ir + 
invariant mass distribution in the data superimposed with 
the background generated according to our procedure. There 
is now a clear agreement with the data distribution, except 
for the K* + peak which the fake pair technique cannot re- 
produce. 

We have checked the background generation procedure 
on samples of A — > pir~ and Kg — > ir + ir~ events, using 
only tracks originating from the primary interaction ver- 
tex in order to artificially increase the background. Fig. [5] 
shows pix~ and 7r + 7r~ invariant mass distributions super- 
imposed with the predicted background estimated using 
our fake pair procedure. From these plots we conclude that 
our background generation procedure provides a realistic 
background estimate in these cases as well. 

It is worth mentioning however that the fake pair tech- 
nique does not take into account contributions from res- 
onances which introduce correlations in M different from 
those generated by making random pairs. 17(1660)— 17(2250) 
decaying into proton and kaon and K* — ► K s ir decays 
with the pion taken as the proton might be important 
sources of distortion of the background shape. With help 
of MC we find a negligible contribution of £'(1660) — 
£(2250) resonances to K s p invariant mass distribution 
while K* — > Kgir decays increase the background by 5- 
10% at small M ~ 1500 MeV/c 2 . Therefore we renor- 
malize the background distribution obtained by fake pair 
procedure by a ratio of two K s p distributions obtained 
with help of MC with and without K* — > K s ir decays. 
Finally the obtained background distribution is normal- 
ized to the data distribution at M > 1650 McV/c 2 . Fig. © 
shows the invariant mass distributions of combinations of 
a positively charged track, assumed to be a proton, and a 
Kg for the data and for the fake pair background, with- 
out using proton identification and with "optimal " proton 
identification. The "signal" interval 1510 < M < 1550 
McV/c 2 is excluded in the data. There is good agreement 
between the shapes of the data and background distribu- 
tions. Polynomial fits of the data excluding the "signal" 
interval 1510 < M < 1550 McV/c 2 is also shown as dashed 
curves. There is a reasonable agreement of the background 
shapes obtained by fake pair procedure and by a polyno- 
mial fit of the data. 
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5 <9 + analysis tools 

5.1 The proton identification strategy 

The + signal is expected to appear as a narrow peak in 
the invariant mass distribution of Kg -proton pairs. K s 
are identified using their U°-like signature (see Sec. I3.1[) . 
To separate protons from 7r + , for each positively charged 
track we build likelihoods under the proton and tt + hy- 
pothesis using the information from DC, TRD, and ECAL 
(see Sec. 3.2), and we take their ratios Cdc, £trd, C-ecal' 



Ctrd , £ecal simultaneously in each interval and find 
those cuts which maximize the signal / 1 ^/background 
ratio. 

We check this procedure on a sample of A — > pir~ events. 
Fig. [7] displays the invariant mass distributions of proton- 
7r~ pairs in both MC and data without proton identifi- 
cation and with "optimal" for the A — > pir~ observation 
proton identification, for —0.6 < xp < —0.3. With "opti- 
mal" proton identification the significance of the A — ► p7r~ 
signal increases in both MC and data samples. 



Cdc{p,L), L— track length 

£trd(p, ctrd), ctrd - energy release in TRD 
£ecal{p, cecal), cecal - energy release in ECAL 
p — track momentum. 

We optimize the cuts for the proton identification like- 
lihood ratios maximizing the sensitivity to the expected 
+ signal. These "optimal" cuts are not necessarily those 
which maximize the purity of the proton sample. 

The best approach for tuning the proton identification 
cuts would be to maximize the sensitivity using a detailed 
Monte Carlo for <9 + production. However, given the poor 
knowledge on the properties of this particle, there is no 
available MC generator describing the production of ex- 
otic baryons. We create, therefore, "fake" + states in the 
NOMAD event generator by using pairs of protons and Kg 
with invariant mass close to the mass of + state. How- 
ever, in this approach the momentum distribution of these 
"fake" + states is determined by the momentum distri- 
bution of protons and Kg from the primary vertex. This 
can result in wrong "optimal" cuts if the true momentum 
distribution of + particles is very different. We try to 
avoid this problem by subdividing the original MC sample 
into several narrow bins of xf and optimizing the cuts for 
each x p interval independently. The xf variable is defined 
as the ratio of the longitudinal projection of the + mo- 
mentum on the hadronic jet momentum to the hadronic 
jet energy in the hadronic center-of-mass frame. The vari- 
able xf is in the range (—1,1) with negative (positive) 
values often called the target (current) fragmentation re- 
gions. 

The procedure of tuning the proton identification cuts 
is then as follows: 

— We build "fake" <9 + states by taking AT^-proton pairs 
with 1510 < M < 1550 MeV/c 2 . Assuming no 0+ po- 
larization, a flat distribution of cos*?*, where 9* is the 
angle between the proton momentum in the + rest 
frame and the + momentum in the laboratory. We 
reweight the cos 9* distribution so obtained to make it 
flat. This is our MC "signal". 

— Any other combination of a K s and a positive track 
not identified as a proton, but with an assigned proton 
mass, is taken as the MC background if its invariant 
mass M falls in the same mass interval. 

— We split the "fake" + states into several intervals of 
positive track momentum. We vary the cuts on Cdc, 





N(pK%) 


purity (in %) 


all 


"signal" 


all 


"signal" 


no ID 


53463 


1856 


23 


16.4 


"optimal" ID 


40561 


1090 


27.8 


22.1 



Table 2. Numbers of pK% pairs and purity of proton sam- 
ples in the data for two subsets of events: without proton 
identification and with "optimal" proton identification. These 
numbers are shown for all entries and for "signal" region: 
1510 < M < 1550 MeV/c 2 . 



In Tab. [2] we show numbers of pKg pairs and purity 
of proton samples in the data for two subsets of events: 
without proton identification and with "optimal" proton 
identification. These numbers are shown for all entries and 
for "signal" region (1510 < M < 1550 MeV/c 2 ). 



5.2 The pKg mass resolution 

The expected mass resolution of the pK s pair is estimated 
as follows. 

— For MC events we calculate the invariant masses of 
the generated and reconstructed pKg pairs, and we 
fit the distribution of the difference between the two 
values by a Gaussian whose width is taken as the mass 
resolution (method "A"). 

— Using the measured momenta of the proton (pi) and 
of the Kg (P2), the angle 9 between pi and P2, and the 
associated errors cr(pi) and cr(p2) we find (neglecting 
errors in cos 9): 



( |rPi-P2 cos9 



—P2 - Pi COS9 ) <T 2 (p 2 )- 



(3) 



This method, "B", can be applied to both MC and data 
events. 

Fig. [8] displays the expected mass resolution of pKg 
pairs as a function of their reconstructed invariant mass, 
as obtained using method "A" (MC only), or method "B" 
(for both MC and data) . The results agree well with each 
other and predict a resolution of about 8.8 MeV/c 2 at the 
0+ mass (1530 MeV/c 2 ). 
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5.3 The statistical analysis 

An estimation of the signal significance in the data is per- 
formed as follows: 

1. A possible difference in the proton cos 9* distribution 
for the signal and background is exploited to improve 
the signal sensitivity. We take all i^-proton pairs with 
1510 < M < 1550 MeV/c 2 , and wc split them into 
10 intervals with similar statistics: five mass intervals 
with cos#* in the interval [—1, —0.5), and another five 
mass intervals with cos#* in the interval [—0.5, 1]. The 
total mass interval (1510 < M < 1550 MeV/c 2 ) covers 
well the expected + mass. The mass bin width, 10 
MeV/c 2 , is comparable to the expected invariant mass 
resolution of K s -proton pairs. 

2. We compute two likelihoods: 

In Lb = 2J + n i ' m M 

1=1,10 , . 

(4) 

lnL B +s= (— h - Si + Hi ■ In (bi + Si)) 

i=l,10 

where hi, Si, rii are the number of predicted back- 
ground and signal events, and observed data events 
in the i-th bin. 

3. We compute the signal statistical significance as: 

S L = ^2{\nL B+s -\nL B ) (5) 

4. We find the resonance mass position M and Brcit- 
Wigner width r and the number of signal events N a 
which maximize Sl- 

For the background we use the procedure decribed in Sec.[U 
The signal is modeled by a Breit-Wigner distorted by a 
Gaussian resolution with a = 8.8 MeV/c 2 . This algorithm 
was checked on several generated distributions containing 
a Breit-Wigner signal of width r distorted by a Gaus- 
sian resolution of width a and superimposed on a fluc- 
tuating background. We considered three cases, a <C F , 
a = r, a ^> r, and found that in all cases the proce- 
dure of maximizing Sl correctly determined the number 
of signal events and r (with r around zero for the case 

tr»r). 

5.4 Opening the box 

We split the data into five xp intervals: 
[-1,-0.6), (-0.6,-0.3), (-0.3,0), (0,0.4), (0.4, 1]. In each 
interval we optimize the proton identification cuts as de- 
scribed in Sec. 13.21 and estimate a possible signal in the 
region 1510 < M < 1550 MeV/c 2 as described in Sec.lOl 
Figs. [9]fT3l display the results. From these plots we con- 
clude that we observe no evidence for the + state in any 
x f interval. In Fig. [H] we display the invariant mass dis- 
tributions of combinations of a positively charged track 
(assumed to be a proton) and a Kg for the two cases of 
no proton identification and optimum proton identifica- 
tion, for —l<xp<l. Table [3] summarizes the results 



and provides also the upper limits at 90% confidence level 
(CL) on the number of + s candidates (-/V" p ) and on the 
production rate R up for both cases. The calculation of the 
upper limits for the production rate include corrections 
for inefficiencies, including the lack of detection of 
mesons, and take into account the K s — ► tt + it~ branch- 
ing ratio. The results arc presented for each bin of xp, 
and also integrated over xf. Fig. [TB] displays the sensi- 
tivity and upper limits (90% CL) for the + production 
rate as a function ot xf- The upper limits are given as five 
curves, each corresponding to a fixed + mass, obtained 
by varying both the number of signal events and the + 
width to maximize Sl as outlined in Sec. 15.31 

We also measure the xf distribution of a potential + 
state as follows. We build the xf distributions in two side- 
bands, 1460 < M < 1500 MeV/c 2 and 1580 < M < 1600 
MeV/c 2 . Wc then normalize the average of these two dis- 
tributions to the expected number of background events 
in the "signal" region (1510 < M < 1550 MeV/c 2 ), and 
subtract it from the xp distribution of the data in the 
"signal" region. The result can be considered as the xf 
distribution of the signal, and could shed a light on the + 
production mechanism. Fig. 1 151 displays the result with no 
proton identification and with optimal proton identifica- 
tion. We observe no statistically significant accumulation 
of events at any xf value. 

6 Conclusions 

We have performed a blind search for the + exotic baryon 
in the 0+ -¥ p + K s decay mode in the NOMAD v^N 
data. We have built a robust background estimation pro- 
cedure which has been tested against various known cases 
like A — ► p7r~, Kg — > 7r + 7r~ and K* — > K s ir. In all cases 
good agreement between data and estimated background 
has been found. Good agreement has also been found be- 
tween the invariant mass (M) distribution of Af^-proton 
pairs in the data and the estimated background in the 
whole mass region excluding the + signal region. We 
have developed proton identification tools based on the 
discrimination power of three sub-detectors, and we have 
tuned the proton identification criteria by maximizing the 
sensitivity to the expected signal in five Xf intervals in- 
dependently. We have checked this approach for A — ► p7r~ 
and found that this procedure indeed maximizes the sig- 
nal significance in both MC and data. Finally, we have 
"opened the box", i.e. examined the + signal in the 
data and found good agreement between the data and the 
background for the whole M region, including the "signal" 
region, in each xf interval. We observe no evidence, there- 
fore, for any + signal in the + — > p + K s decay channel 
in the NOMAD v^N data. We give an upper limit at 90% 
CL on + production rate of 2.13 ■ 10 -3 events per neu- 
trino interaction at M = 1530 MeV/c 2 after integrating 
over xf- 

It is interesting to compare this result with the re- 
cent analysis of old bubble chamber neutrino experiments 
which provide an estimation of the + production rate 
as large as ~ 10 -3 events per neutrino interaction [10] . 
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if interval 


[-1,-0.6) 


(-0.6,-0.3) 


(-0.3,0) 


(0,0.4) 


(0.4, 1] 


all 


no ID 














N s (fit) 


18 


26 


35 


30 


65 


77 


SL 


1.96 


1.49 


1.01 


1.18 


2.61 


1.82 




41 


61 


88 


81 


101 


161 


R up 


3.84 


2.18 


1.74 


1.37 


0.83 


4.36 


optimal ID 














N s (fit) 


12 


29 


-26 


-34 


24 


-33 


SL 


1.38 


1.72 


1.35 


1.85 


1.25 


0.97 


AT 


28 


68 


39 


36 


52 


67 


R up 


2.80 


2.60 


0.84 


0.79 


1.00 


2.13 



Table 3. Upper limits (90% CL) on the number of + candidates (N™ p ) and on the <9 + production rate (R up , in units of 
events per 10 3 interactions) for the case of no proton identification and with optimal proton identification. 



As shown in FigfTol for a large fraction of the xp range, 
except in the region xp — > — 1, such a value is excluded. 
Unfortunately, ref. [lOj does not provide information on 
the xf region in which a + signal was observed. Fur- 
thermore, in ref. |10| we find no information that the back- 
ground estimation procedure took into account the effects 
mentioned in Sec. [4J which can result in an underestima- 
tion of the background and thus in an overestimation of 
both the signal significance and the production rate. 

Preliminary NOMAD results from searches for the ex- 
otic + baryon reported earlier [3B], quoting a hint for 
a signal with a statistical significance of 4.3 a, suffered 
from an incorrect background estimation, which did not 
take into account the effects mentioned in Sec. 21 The re- 
sults reported in |36| are obtained on a smaller sample 
of the NOMAD data. The positives of that sample were 
subjected for a cleaner proton identification which yielded 
an increase of the purity of the protons sample from 23% 
to 51.5% with about factor six lost of the statistics. The 
difference in shapes of K^p invariant mass distributions 
reported in [3^1 and in Figs. [5] is due to an additional 
requirement imposed in [36j on energy of protons to be 
larger than that of Kg. 
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Fig. 3. The mean value of the invariant mass of pn~ pairs 
from identified A's as a function of the reconstructed proton 
momentum with no momentum correction (full circles) and 
with momentum correction (stars). 




Fig. 4. Kgir invariant mass distribution in the data, superimposed with the background generated without hadronic jet 
rotation (left), and with hadronic jet rotation (right). 
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Fig. 6. Invariant mass distribution of pairs of one positively charged track (assumed to be a proton) and a Kg for the data 
(points with error bars) and for the fake pair background (shadowed area). (Left) Proton identification has not been used. 
(Right) "Optimal" proton identification. Data in the "signal" region (1510 < M < 1550 MeV/c 2 ) are not shown. Dashed curve 
is fit of the data by a polynomial excluding the "signal" region. 
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Fig. 7. Invariant mass distributions of proton-7r~ pairs for —0.6 < xf < —0.3 interval. (Left, up) MC, no proton identification; 
(Right, up) MC, with "optimal" proton identification; (Left, down) data, no proton identification; (Right, down) data, with 
"optimal" proton identification. The curves represent the predicted background and the amount of A — > pn~ signal maximizing 
the signal / ^background ratio (see Sec. 5.1). 
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Fig. 8. Expected invariant mass resolution of proton-JsTg pairs 
as a function of the invariant mass for method "A" (MC only), 
and for method "B" (MC and data). See text for details. 
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Fig. 9. Invariant mass distributions of pairs of a positively charged track (assumed to be a proton) and a K% in the data, for 
— 1 < xf < —0.6. (Left) No proton identification; (Right) Optimized proton identification. The curves represent the predicted 
background and the amount of & + signal maximizing Sl (see Sec. 5.3). 
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Fig. 10. Invariant mass distributions of pairs of a positively charged track (assumed to be a proton) and a K% in the data, for 
—0.6 < if < —0.3. (Left) No proton identification; (Right) Optimized proton identification. The curves represent the predicted 
background and the amount of + signal maximizing Sl (see Sec. 5.3). 
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Fig. 11. Invariant mass distributions of pairs of a positively charged track (assumed to be a proton) and a K% in the data, 
for —0.3 < xf < 0. (Left) No proton identification; (Right) Optimized proton identification. The curves represent the predicted 
background and the amount of + signal maximizing Sl (see Sec. 5.3). 
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Fig. 12. Invariant mass distributions of pairs of a positively charged track (assumed to be a proton) and a K% in the data, 
for < xf < 0.4. (Left) No proton identification; (Right) Optimized proton identification. The curves represent the predicted 
background and the amount of + signal maximizing Sl (see Sec. 5.3). 
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Fig. 13. Invariant mass distributions of pairs of a positively charged track (assumed to be a proton) and a K% in the data, 
for 0.4 < if < 1. (Left) No proton identification; (Right) Optimized proton identification. The curves represent the predicted 
background and the amount of + signal maximizing Sl (see Sec. 5.3). 
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Fig. 14. Invariant mass distributions of pairs of a positively charged track (assumed to be a proton) and a K% in the data, 
for — 1 < xf < 1- (Left) No proton identification; (Right) Optimized proton identification. The curves represent the predicted 
background and the amount of + signal maximizing Sl (see Sec. 5.3). 
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Fig. 15. xf distribution of a potential + signal in the data. (Left) No proton identification; (Right) Optimized proton 
identification. 
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Fig. 16. Sensitivity and upper limits at 90% CL for + pro- 
duction rates as a function of xf, for + masses of 1510, 1530, 
1550 MeV/c 2 . 



